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Abstract: The spin-torsion theory is a gauge theory approach to gravity that expands upon Einstein’s 
general relativity (GR) by incorporating the spin of microparticles. In this study, we further develop 
the spin-torsion theory to examine spherically symmetric and static gravitational systems that involve 
free-falling macroscopic particles. We posit that the quantum spin of macroscopic matter becomes 
noteworthy at cosmic scales. We further assume that the Dirac spinor and Dirac equation adequately 
capture all essential physical characteristics of the particles and their associated processes. A crucial 
aspect of our approach involves substituting the constant mass in the Dirac equation with a scale 
function, allowing us to establish a connection between quantum effects and the scale of gravitational 
systems. This mechanism ensures that the quantum effect of macroscopic matter is scale-dependent 
and diminishes locally, a phenomenon not observed in microparticles. For any given matter density 
distribution, our theory predicts an additional quantum term, the quantum potential energy (QPE), 
within the mass expression. The QPE induces time dilation and distance contraction, and thus mimics 
a gravitational well. When applied to cosmology, our theory yields a static cosmological model. 
The QPE serves as a counterpart to the cosmological constant introduced by Einstein to balance 
gravity in his static cosmological model. The QPE also offers a plausible explanation for the origin 
of Hubble redshift (traditionally attributed to the universe’s expansion). The predicted luminosity 
distance-redshift relation aligns remarkably well with SNe Ia data from the cosmological sample of 
SNe Ia. In the context of galaxies, the QPE functions as the equivalent of dark matter. The predicted 
circular velocities align well with rotation curve data from the SPARC (Spitzer Photometry and 
Accurate Rotation Curves database) sample. Importantly, our conclusions in this paper are reached 
through a conventional approach, with the sole assumption of the quantum effects of macroscopic 
matter at large scales, without the need for additional modifications or assumptions. 


Keywords: alternative theories of gravity; dark matter; dark energy; galactic rotation curves 


1. Introduction 


The primary challenge facing the current theory of gravity, the Einstein-Newton 
theory, lies in the unresolved mysteries of dark matter and dark energy. To date, there has 
been a lack of direct observational evidence confirming the existence of these enigmatic 
components. Extensive efforts are being made to address this topic. On the theoretical front, 
much attention is directed towards modifying or extending the Einstein—Newton theory of 
gravity to align with astronomical observations. For instance, one can enhance the standard 
Lagrangian in general relativity by incorporating higher-order curvature corrections [1-6], 
or formulate non-linear Lagrangians [7,8]. Other relevant examples include modified 
Newtonian dynamics (MOND) [9,10] and its relativistic version [11], as well as conformal 
gravity [12,13]. 

In this paper, we explore the concept of quantum spacetime at cosmic scales as a 
potential alternative to dark matter and dark energy. This research represents a significant 


Universe 2024, 10, 333. https:/ /doi.org/10.3390/universe10080333 


https: //www.mdpi.com/journal/universe 


Universe 2024, 10, 333 


2 of 38 


extension of our previous study [14]. Rather than relying on previous assumptions, we 
explicitly propose scale-dependent quantum properties of spacetime. This proposition is 
motivated by our interest in potentially moving away from the notion of a preferred abso- 
lute spacetime (PAS) for the universe, as implied by general relativity (GR) and illustrated 
in the standard ACDM cosmology. It is believed that the PAS is theoretically characterized 
by the Friedmann—Lemaitre-Robertson—Walker (FLRW) metric and empirically supported 
by the cosmic microwave background (CMB). One aim of this study is to demonstrate 
that the notion of a PAS is not inherently valid when considering the quantum effects of 
macroscopic matter at large scales. 

Before going further, it is imperative to clarify our actual motivation in order to prevent 
any possible misunderstandings. For ease of reference, we categorize the universe into 
three distinct classes based on spatial scales: the cosmic scale (encompassing galaxies to 
the entire universe), the macroscopic scale (ranging from everyday life to the size of the 
solar system), and the microscopic scale (where quantum mechanics becomes significant, 
as is well-known). According to the mainstream view, the quantum effect for objects at the 
macroscopic and cosmic scales can be safely disregarded, due to their large mass. This is 
known as the classical limit. We will show, however, that quantum effects are dependent 
on the scale of the gravitational system considered and that the “classical limit” is only 
reached at macroscopic scale. In other words, quantum effects are significant not only at 
microscopic scale but also at cosmic scale, and can only be ignored at macroscopic scale. 
Intuitively, when we study cosmology, any observers who perceive distant galaxies in the 
universe, would expect them to exhibit quantum behavior resembling that of microscopic 
particles, due to their distant location. As for galactic dynamics, quantum effects should 
also be taken into account for galaxies, although in the central part, they can be ignored 
due to the mechanism presented in this study. 

It is helpful to discuss the reasons behind generalizing the quantum effects observed in 
microscopic particles to macroscopic ones, as well as why the strength of these generalized 
quantum effects appears to increase with the spatial scale, based on our previous paper [14]. 
It is believed that gravity turns disorder into order, and order is fundamental to space, time, 
and spacetime—or, as we might say, spacetime is nothing but the continuous ordering of 
events. Thus, without gravity, there would be no spacetime. In any local inertial frame, 
gravity is still present, but the net gravitational force is canceled out by the inertial force, 
leaving order or flat spacetime behind. In special relativity, rigid rods are used to create 
coordinate lattices, ensuring that every event has a position at any given time. According to 
GR, the global structure of curved spacetime can still be described using arbitrarily curved 
rigid rods (depending on the matter distribution) of arbitrary length. However, quantum 
mechanics reveals that microscopic particles can escape the order depicted by spacetime 
and thus elude the control of gravity. The quantum effects of electrons, for instance, can be 
attributed to their quantum randomness (also known as intrinsic randomness), which is 
markedly different from classical randomness (or apparent randomness), exemplified by 
Brownian motion. A free electron deviates from a straight line or a curved geodesic, which 
implies that gravity is not strong enough to bring it into order. This explains quantum 
effects from a geometric point of view. Gravity does bring macroscopic matter into order 
in the sense that any free macroscopic object will move along a straight line or a curved 
geodesic. This fact is fundamental to both Newtonian theory and GR. However, if we 
assume that all matter, irrespective of its mass, has two opposite properties—gravity 
and quantum randomness—then it is possible that gravity is not strong enough to bring 
any distant macroscopic object into order, just as occurs with microscopic particles. If 
our assumption is true, then general relativity (GR) is valid only approximately within a 
sufficiently small neighborhood of any point in the spacetime manifold. In this case, the 
uncertainty in matter distribution, originating from quantum randomness, will accumulate 
with spatial distance. Or, put another way, the strength of quantum effects appears to 
increase with spatial scale. On the other hand, since the spin-induced torsion arising from 
macroscopic matter vanishes locally, the geodesics of test particles are precisely defined. 
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This serves as the physical foundation for physicists to construct coordinate lattices. As 
a result, in general, spacetime is not only curved but also flexible. It turns out that the 
preferred absolute spacetime (PAS), which is determined by all matter in the universe, loses 
its meaning when scale-dependent quantum effects come into play. Any new fundamental 
assumptions about physical laws should be falsifiable. The rest of this study represents a 
first attempt along this line. 

We will illustrate the integration of our assumption with the Einstein equations. It is 
natural for us to generalize the Dirac theory to describe the macroscopic matter content of 
gravitational systems with cosmic scale. 

There are two approaches to interpreting quantum mechanics. The first is the standard 
view, which asserts that microscopic particles cannot have continuous trajectories and that 
non-commuting observables (such as position and momentum) must satisfy the uncertainty 
principle. The second approach is known as the causal interpretation, which posits that each 
microscopic particle has a continuous trajectory and replaces the uncertainty principle with 
the concept of quantum potential [15,16]. Both approaches yield equivalent predictions 
for experimental results. However, the causal interpretation, particularly in its geometric 
algebra (GA) version, provides profound insights into the quantum nature of spacetime and 
matter [17-20]. In this paper, similarly to in our previous work, we utilize spacetime algebra 
(STA) [20-22] as our mathematical language. STA is constructed from a Minkowskian vector 
space and provides a straightforward geometric understanding of Dirac theory. It provides 
a natural link to classical mechanics, and for expressing Dirac theory with observables, it 
offers enhanced computational efficiency and capabilities [18] when compared to the tensor 
analysis method [23]. The fundamentals of STA are outlined in Appendix A. In the STA 
version of Dirac theory, known as real Dirac theory, Hestenes identified the complex number 
present in the matrix version of Dirac theory as the spin plane [17]. Furthermore, since the 
complex number is always accompanied by the Planck constant fh, a rigorous derivation of 
the Schrodinger theory from the Pauli or Dirac theory implies that the Schrodinger equation 
describes an electron in an eigenstate of spin, rather than, as commonly believed, an electron 
without spin [24]. Correspondingly, there are two gauge theories of gravity concerning our 
research work: one is the tensor analysis approach [25], and the other is the gauge theory 
gravity (GTG) developed by the Cambridge group using STA [26]. These two theories are 
nearly equivalent, but the latter is conceptually clearer and technically more powerful for 
comprehending and calculating the challenges encountered in quantum mechanics and 
gravity. Therefore, we choose to adopt GTG in this paper. The fundamentals of GTG and 
its applications to spin-5 particles are summarized in Appendix C. 

It is now widely acknowledged that the quantum random motion of spin-4 particles 
can be fully described by their spin. In the presence of gravity, the spin gives rise to 
the torsion of spacetime [21,25—29]. In particular, the effects of spin-torsion in GTG were 
investigated in Ref. [27]. GTG is nearly equivalent to Einstein—Cartan gravity [25]. The 
stress-energy tensor derived from the Dirac theory contains an asymmetric component, 
representing the contribution of quantum spin. As a consequence, the metric of spacetime 
as predicted by the generalized Einstein equations includes a component that accounts 
for the torsion of spacetime. As will be shown later, the torsion term appearing in the 
metric represents a modification to GR. To date, the problem of finding solutions for a 
Dirac field coupled to gravity in a self-consistent manner has primarily been considered 
for microscopic particles. To our knowledge, there is only one work that explored the 
massive, non-ghost cosmological solutions for the Dirac field coupled self-consistently to 
gravity [30], which is close to our present study. As expected, the authors applied their 
methods to the very early universe, with the Dirac field describing massive yet microscopic 
particles. This differs from our present study, in which we assume that the Dirac field 
describes macroscopic particles within a static universe. 

For massive fermions, such as electrons and neutrons, the mass will always appear 
in the phase factor of the solutions of wave equations. When gravity matters, this mass 
dependence remains. Consequently, it is generally believed that the gravity effect is not 
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purely geometric at the quantum level. Hence, the immediate question is how can we 
reconcile, if possible, the contradiction between the quantum effect and the requirement of a 
geometric description of gravity? For microparticles, reconciling this contradiction is indeed 
impossible. However, in this study, we propose that the scale-dependent quantum effect 
on macroscopic matter only becomes significant at cosmic scales. Locally, at a macroscopic 
scale, the large mass m of each particle, or equivalently, the small value of the number 
density p of the fluid, guarantees what is known as the classical limit, allowing for a 
geometric description of the gravitational effect on macroscopic matter. The remaining 
question is what mechanism ensures that quantum effects are significant at cosmic scales 
while being negligible at macroscopic scales? 

The answer to this question is not intricate, but rather subtle. The specifics will be 
provided later, but for now, we will outline the main concepts behind the mechanism. 
Our analysis commences with an examination of the Dirac theory within various inertial 
frames [18]. Generally, the Lorentz—invariant Dirac spinor is defined in spacetime as 


p(x) = p'/7eB)/2 R(x). (1) 


where p(x) is a scalar representing the proper probability density and R(x) is a rotor 
(Lorentz rotation) satisfying RR = 1. In Dirac theory, the parameter A(x) is intriguing. 
It is noteworthy that the states of a plane-wave particle have $ = 0, while those of an 
antiparticle have 6 = 7. In this paper, we aim to generalize the spinor y(x) to describe 
free-falling macroscopic, non-relativistic particles with identical mass m. Therefore, it is 
logical to set 6 = 0 throughout. We thus adopt [14] 


p(x) = p(x)? R(x). (2) 


For our purpose, we can interpret pọ = py) as the proper number density of a fluid, then we 
refer to Pm = myi as the corresponding proper mass density. As illustrated in Appendix B, 
the rotor R can be used to transform a fixed frame {7} into a new frame {e, = Ry,R}, 
and we identify v = eg as the proper velocity associated with the expected history x(T) of 
a particle and thus the current velocity of the fluid. Naturally, the corresponding stress- 
energy tensor of the fluid can be written as 


T(a) = Pma- vv. (3) 


It is important to point out that this stress-energy tensor, which originates from the Dirac 
spinor given by Equation (2) for a single electron, definitely describes a classical pressure- 
free ideal fluid without incorporating quantum spin. Quantum spin is only included if 
the spinor satisfies the Dirac equation, as will be discussed soon. In this sense, we can say 
that the spinor given in Equation (2) captures all aspects of macroscopic particles, when 
considered with or without the Dirac equation. 

The Dirac equation is 


AV piy3 = my. (4) 


where fi is the Planck constant. This equation describes a single spin-5 free particle with a 
fixed mass m and a probability density p. It is convenient to define a spin density trivector 
as (note that this is denoted with S3 in Appendix B) 


ae ee oe 
S = J PtysP = FPRIYR. (5) 


We show in Appendix B that, from the Dirac Equation (4), the general stress-energy tensor 
for the spinor field is 


T(a) = pma -vv + [a - V (Sv)] -v — (aA V)- (Svv. (6) 


Universe 2024, 10, 333 


5 of 38 


When S = 0 (this is the classical limit), this equation becomes the stress—energy tensor of a 
pressureless ideal fluid without spin, as shown in Equation (3). Clearly, if we treat each 
free-falling macroscopic particle with a fixed mass m as a single spin-} free particle, we 
have the opportunity to include its quantum behavior under certain circumstances. This is 
the main aim of this study. 

We are now in a position to discuss how the spinor approach can unify quantum 
mechanics and classical mechanics. We believe this will be helpful to understand the 
physical significance of our proposal in this study. The general form of the Dirac spinor, 
as defined in Equation (1), and the corresponding Dirac equation are translations of the 
traditional matrix version of Dirac theory, expressed in terms of STA [18]. They are equiva- 
lent in practical applications. The advantages of the STA version of Dirac theory include 
providing an explicit geometric and causal interpretation of the theory, which makes it very 
convenient when incorporating with general relativity (GR). One of the unique features of 
the Dirac equation is that it allows for solutions with both positive and negative energy 
states. The negative energy solutions in the Dirac equation lead to the concept of the Dirac 
sea, a theoretical model that was used to explain the existence of antiparticles. Because of 
the presence of negative energy states, when dealing with Dirac spinors, one encounters 
different types of physical quantities called densities, such as scalar densities, which are 
Lorentz-invariant. Scalar densities can interact with scalar fields, and the Higgs field is one 
such example, which is responsible for giving mass to particles through their interaction 
with it. When we generalize the Dirac theory to describe classical macroscopic particles, 
however, we restrict ourselves to the positive energy states. This is adequate because, in the 
case of the non-relativistic energy of macroscopic particles, there are no particle—antiparticle 
creation processes involved. That is, we treat each macroscopic particle as a single Dirac 
fermion. As such, we recognize the proper number density p = py and the proper mass 
density Pm = mip (as measured by observers who are comoving with the current velocity 
v = Ryok) as that for the ensemble of a single macroscopic particle, described by the Dirac 
spinor give in Equation (2). As mentioned, classical mechanics are recovered when S = 0. 

When applying these results to gravitational systems, we simply need to replace V 
with the covariant derive symbol D, while all other quantities remain gauge-invariant. 
As such, the classical limit (also known as the short wave approximation) refers to the fact 
that the contribution of the spin to the stress—energy is negligible. Clearly, the magnitude 
of the spin |S| = Ëp is determined by the number density p. Therefore, an alternative 
interpretation of the classical limit is that the value of the number density, p, is ignorably 
low, whereas a significant quantum effect arises when p is sufficiently high. For a fixed 
mass density pm(x), the spin |S| can be determined via p(x) = Pm(x)/m. This freedom 
strongly indicates that, in order to incorporate the scale-dependent quantum effect at scale, 
we can replace the mass m in the Dirac equation with a mass function that depends on 
the spatial scale A of the system, denoted as m(A). The mass function m(A) naturally 
satisfies the condition that, when A — 0, then S — 0, corresponding to the macroscopic 
scale. Conversely, when A — œ, then S becomes a A-dependent quantity related to the 
constant density of the universe. The intermediate functional form of m(A), which plays a 
crucial role in the study of galaxies, can be determined through observations. Consequently, 
we can employ the Dirac equation in the presence of gravity without any modifications, 
as there are no derivatives of m(A) involved in our calculations. The properties of the mass 
function m(A) only need to be discussed in the final results. 

It is important to explicitly acknowledge that when applying the Dirac equation to 
gravitational systems at cosmic scales, we shall show in Section 3 that the anisotropy of 
the spin density adheres to the symmetric properties of mass density m(x). This fact only 
applies to macroscopic matter according to our assumption, where the gravitational effect 
is geometric throughout the spacetime manifold and, furthermore, the proper mass density 
and the spin density have already been averaged and must satisfy the condition that the 
quantum effect vanishes locally. In contrast, when considering electrons, the anisotropy 
of spin must be taken into account in all cases, as the gravitational effect is not purely 
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geometric. Thus, the quantum effect can never disappear locally, as is commonly under- 
stood [30-33]. 

The subsequent sections present an examination of radially symmetric and static 
gravitational systems based on our new postulation in Section 2, the applications of our 
findings to cosmology and galaxies in Section 3, and summarized conclusions and discus- 
sions in Section 4. We employ natural units (G = A = c = 1) throughout, except where 
stated otherwise. 


2. Radially Symmetric Gravitational Systems 


We want to investigate the quantum nature of spacetime at cosmic scales, with the main 
subjects of concern being galactic dynamics and cosmology. Traditionally, a pressureless 
ideal fluid has been a good model for both galactic dynamics and matter-dominated 
era cosmology. In both scenarios, matter experiences free-fall in a gravitational field 
and is characterized by the density distribution p(x) and velocity field v(x), where x 
represents the position vector in spacetime. This model is usually known as collapsing dust. 
In gravitational systems, it is widely recognized that a static state can only be sustained 
when there is a substantial gravitational potential well, as observed in galaxies. However, 
in the realm of cosmology, the homogeneous and isotropic distribution of matter across the 
entire universe renders static solutions non-existent. 

Nevertheless, the introduction of quantum effects at cosmic scales fundamentally 
alters the situation. 

The Dirac theory for radially symmetric gravitational systems can be investigated 
using gauge theory gravity (GTG). The latter is constructed such that the gravitational 
effects are described by a pair of gauge fields, h(a) = h(a,x) and w(a) = w(a,x), defined 
over a flat Minkowski background spacetime [26], where x is the STA position vector and is 
usually suppressed for short. Luckily, the majority of the necessary results for our present 
study have already been derived in prior work [27]. These results are valid for microscopic 
particles; what we need to do is to generalize these to macroscopic matter. 

Let us consider a radially symmetric and static gravitational system composed of 
free-falling particles with the identical mass m. We make the assumption that the Dirac 
spinor y(r) = p(r)!/2R(r) defined in (2) can fully capture all the essential physical aspects 
of the system if it satisfies the Dirac equation 


Dyiys = my. (7) 


We interpret p = py as the proper number density and pm = mi as the proper mass 
density. We adopt S = 5piy3) = }pRiy3R defined in (5) as the spin density trivector for 
the gravitational system. So there are four variables, namely m, p, Om, and S, that can be 
used to describe the gravitational system. However, only two of them are independent, 
as they must satisfy the conditions p = pm/m and |S| = 4p. In this paper, we eliminate 
m and retain the other three variables. Among them, the relationship between pm and p 
can be further determined through observations by assuming a specific form for the mass 
function m(A). Once this has been carried out, we are left with only one variable, which 
is the proper mass density pm. Therefore, if we are given the proper mass density pm of a 
gravitational system, we can predict the entire set of results, including the quantum effects. 
Now we define a set of spherical coordinates. From the position vector of the 

flat spacetime 
x = tyo +r sin 0(cos py; + sin y2) +r cos 073, (8) 


we obtain the basis vectors, as follows: 
et = OX =Y0, 


er = 0,-x =sin 0 (cos yı + sin y2) + cos O73, 
eg = 0px =r cos O (cos Py; + sin y2) — r sin 673, 


(9) 


ep = Ogx =rsinO(—sin py + cos f72). 
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Since eg and eg are not unit, we define 
6=e9/r, $= eg/(rsin6). (10) 


With these unit vectors, we further define the unit bivectors (relative basis vectors for 
et = Yo) 


Oy =Cret, 
0@ =6e,, (11) 
op =e. 

These bivectors satisfy 

ToT p = CeO =i. (12) 


As in our previous paper [14], we initially attempted to analyze static systems. The set 
of h field that satisfies the spherically symmetric and static matter distribution is assumed 
to take the form [21,26] 

he!) = fie’, he) = ge + ge, 


h(e®) =e, Alet) = e?, 


where f1, g1, and g3 are all functions of r only. We could have tried the form h(e') = fié + fre", 
but it is more reasonable to set fz to zero, which is referred to as the ‘Newtonian gauge’ [26]. 
The subsequent steps of this study can be summarized as follows: 


= DAh(a)=xh(a)-S 
— 


h(e") w(a) 


R(aNb)=Lato(b)—Lyeo(a)+e0(a) xolo) T =3p R(b ^a) \ 
R =0,:R(a) 
SEES, guv =A (u) -KH 
where the metric guy is introduced to compare our results with those predicted by gen- 
eral relativity (GR) and to provide conventional approaches for subsequent applications. 
Notably, w(a) and R(a A^ b) can be decomposed into torsion-free components and torsion 
components [27]. This decomposition allows for the recovery of classical predictions of GR 
when the impact of torsion is insignificant. Conversely, given the extensive research on these 
classical predictions available in the existing literature, we can easily incorporate torsion 
terms into classical results when we deem them to be significant. This allows us to leverage 
decompositions and enhance our understanding of the phenomena under consideration. 
By solving Equation (A154), we can obtain a solution for w(a) [26,27], the result is 


w(a) = w' (a) + sa -S, (14) 


where 


; 
—1 R A a2 
_ (H= )a dð- ($2 )a- doer (15) 


denotes the torsion-free component of w(a), and the new functions G and F are also all 
functions of r only. From the w field, for any bivector B, its strength tensor can be obtained 
directly from Equation (A149) [26,27], the result is 


R(B) =R' (B) + Tee -S)- S 
16) 

1 ( 
— 5K (B -D)-S, 
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where 


R'(B) =a10;B - 0; + (209 + agiog)B - og 
+ (e209 — &3idg)B - op + a607(B A 07) 
+ (#409 — asiog)(B A 09) 
+ (a40p + as5iog)(B A og) 


(17) 


denotes the torsion-free component of R(B), with the understanding that o, \ og = op ^ 
dg = 0, and ay,...,a6 are given by 


w =L- Li + G* — F, 
ht E 
_7 (81 8182 82 
OE E À 
r r r 
wot() +E -(B) 
ws =(—33 +81 —1)/?°. 


From Equations (16) and (17), the Ricci tensor R (a) and Ricci scalar R are given by 
R(a) =R'(a) + 5x2(a-8) -S 
— Ika -(D-S), (19) 


R=R'+ 1282, 


where 
R' (a) =[(a1 + 2&2)a - et + 2asa - e; ler 
+ [223a - e, — (a, + ANa - ey ey (20) 
(a2 +04 + 06)a-00 
(a2 + a4 + a6)a- bp, 
and 


R! = 201 + 4a2 +404 + 2&6, (21) 


denote the torsion-free components of R(a) and R, respectively. 
In order to solve Einstein Equation (A155), we need to know T(a) given by Equation (A157). 
Noting that y = y(r), we find that 


T (a) =(a- Diya) 
- ee | we 

=(a : h(e")dupiysy + sea pirsyy (22) 

=a - (g1e" + gve'){d-pi7s)1 + w(a)- S. 
The terms (0;i73)1 should be derived from the Dirac Equation (7). Nevertheless, as elab- 
orated in our previous paper [14], the use of this stress-energy tensor prototype can be 
misleading when applied to gravitational systems consisting of macroscopic bodies. This 
is because, as S approaches zero, the tensor T(a) defined in (22) also approaches zero. 


This contradicts our expectation that as S tends to zero, T(a) should represent a tensor 
describing a classical pressureless ideal fluid. 
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An appropriate representation of the stress—energy tensor for our specific needs in- 
volves decomposing it into two components: one that characterizes the classical pressure- 
free ideal fluid and another that accounts for quantum effects. By substituting V with D, 
we can rephrase Equation (A113) and express 7 (a) as, 


T (a) = pma- ere, — (a ^A D) - (S-e¢)e¢ + [a - D(S - e¢)] - et, (23) 


where we have replaced v with e; due to the retained gauge freedom to perform arbitrary 
radial boosts in restricting the h function [21,26]. It is important to note that the choice 
v = €: = Yo simply fixes the rotation gauge in such a way that the stress—energy tensor 
takes on the simplest form; there is no other additional physical content. Furthermore, 
setting v = e; does not imply that the fluid is at rest or that the observers are comoving 
with the fluid. In such a setting, all rotation-gauge freedom has been completely removed, 
as it should be, before one can derive a complete set of physical equations. Note also that 
R(B) deals directly with physically measurable quantities, whereas the algebraic structure 
of the h-function is of little direct physical significance. Therefore, the rotation gauge has 
been fixed by imposing a suitable form for R(B), rather than restricting the form of h(a), 
as discussed in Refs. [21,26]. Consequently, the proper mass density pm (r), potential energy 
density po (r), and the corresponding mass (energy) M(r) presented later acquire a similar 
physical meaning as in GR. 
The Einstein equation is 


1 
G(a) = R(a) — zak = KT (a). (24) 
A direct and efficient approach is to solve the Einstein equation separately for two scenarios: 


when a = e; and a = ey. To be specific, in general, our results need to be derived from the 
following equations: 


et: G(ez) =Kez - T (et), 
er: G (er) =Ke; > T (et), 
et: Gler) =Kez- T (er), (25) 
és Gle =Ker -T e: 


However, we demonstrate that, for our intended purpose, solving the first two equations 
listed above is adequate. From (23), we obtain 


er: T (er) = Pm — (ee AD) - (S- et). (26) 
Since (e; A D) - (S-e¢) = (e \ 0g): [a - D(S - e¢)], we need to calculate a - D(S - e;), as follows: 


a-D(S-e) =(a- DS) -e+ S- (a- Dez) 


=[a-h(e")d,S + w(a) x S] -e +S- [w(a) - e]. (27) 
We thus have 
(e Ada) « [a D(S + e)] =(er Ada) - [S - (w(a) +e) 
=et : {[9a A (w(a) - er)] - S} 
=e;- {[—Go;, + «(S - et)]- S} (8) 
=x(e,-S)? = KS? 
Therefore, 
et: Thet) = Pm — KS?. (29) 


Similarly, we have 
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er: T (er) =[e¢- D(S - e¢)] - (e1 Aer) 
=G(S-e,) +(e: ^er) (30) 

=0. 
On the other hand, for G (e+), we find from (19)-(21) that 
1 
G (er) =R(er) = eR 
=203e; — (244 + &6 + 22S )e (31) 
1 1 
+ 5 (er -S)-S— wher: (D-S). 


Substituting (29)-(31) into ez - G (e+) = xer - T (er) and e, - G (e1) = Ke, - T (e1), we obtain 


2&4 +46 = —K(Om— IKS?) (32) 
x3 = 0. 
Interestingly, from (20) we find 
da A R' (a) = 2050; — 2830; = 0, (33) 
which gives 
a5 = 3 = 0. (34) 


By substituting the expressions of as given by (18) into (32) and (34), our solutions to 
the Einstein equations can be summarized as follows: 


$- ia. 3 
2[g13 ($) + (2) - F(&)] - S =- rlon — 58%), 


r r r2 4 
81) , 8182 82 _ 35 
820r( r ) 1 r2 G r 0, ( ) 
8&2; , 8182 81 _ 
g19,( F ) T r2 F F 0. 


The first and third equations in (35) can be combined to give 


Or [r((g2)? — (¢1)? + 1)| = «(om — iss?) r. (36) 


Now, if we define 
r 
M = 5((g2)*— (8: +1), (37) 


we find (remembering x = 872): 


M(r) = 470 f (oni! = FS) ) rar (38) 


The expression of M strongly implies that it can be identified as the mass (total energy) 
of a gravitational system within r at any given time t. Remarkably, the mass M contains 
the quantum effects, —}xS?(r), as intended. Naturally, when S = 0, the conventional 
expression for M in general relativity is regained. 

As mentioned earlier, there are three functions in the } field, namely f1, g1, and g7, that 
need to be determined by solving Einstein equations. So far, we have obtained expressions 
for g1 and gz, and they are linked to the observables pm and S. Thus, the determination of 
fi remains pending. To achieve this, we substitute a = e! into the torsion Equation (A154), 
yielding the subsequent equations 
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DAh(e') =xh(e') -S, 
da A [a-h(e")dy fie’ + frw(a) e] =Kfier-S, 


39 
S19 fi =— Gf, a 
f ee Edr 
On the other hand, from the last two equations of (35), it is easy to derive that 
G = 9;gı, F = drga. (40) 
So, we immediately obtain that 
fi =1/81. (41) 
Now, we turn to the metric tensor Quy = gu ` 8v = h`! (eu) -h-1(e,). From (13) and (41), 
we obtain i 
h! (ei) = 4e; — 26r, h`! (e = —e,, 
a (er) 81et — §2er, N (er) a (42) 


h`! (eg) = eo, h! (ep) = eg. 


So, the metric tensor depends only on g4 and gz. However, the (g3)? included in (37) served 
as the kinematic energy in gravitational systems [26]. To see this, we consider a radially 
free-falling particle with v = e;. We have 


>. dx, . 1 
y= I ter + te, = h(v) = h (et) a + Bp. (43) 


Clearly, in this case, g2 = t represents the radial velocity of the particle. In general, we can 
understand the physical significance of g2 by rewriting (37) as [26] 


se- = 5 ((g)?-1), (44) 


which is a Bernoulli equation for zero pressure and total non-relativistic energy 5((g1)? — 1). 
But for non-relativistic matter, the contribution of kinematic energy to gravity can be 
safely disregarded. So, from (37) we have 


= (1 = ul ae (45) 


r 


This paper specifically concentrates on non-relativistic matter within static gravitational 
systems, enabling us to express the metric as follows: 


-1 
i= (1 MO) ae (1 ee dr? — r2(d0? + sin? Odo”). (46) 


Remarkably, we have obtained a metric that bears a striking resemblance to the familiar 
form in general relativity. However, it is important to note that this metric was derived 
from Dirac theory, and the mass M incorporates the quantum effects arising from macro- 
scopic matter. 

In this paper, the term —}xS?(r) in the mass-energy expression (38) is referred to as 
the “quantum potential energy”. This concept was first introduced by Bohm [34] in 1952 
and later by DeWitt [35] in the same year, with both authors commonly referring to it as 
the “quantum potential” [29]. Similarly to how traditional mass—energy curves spacetime, 
quantum potential energy distorts spacetime. This phenomenon is known as “spin-torsion” 
theory [25,26]. By explicitly decomposing the stress—energy tensor 7 (a) into a spin-free 
part and a spin part, as illustrated in (23), the gravitational strength, represented by the 
Riemann tensor R(B), can be explicitly decomposed into a torsion-free component and a 
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torsion component, as shown in (16). This results in a two-component metric as outlined 
in (46). Naturally, when S = 0, the metric reduces to one that describes a spherically-static 
gravitational system, as suggested by GR. 

It is crucial to grasp the physical significance of the quantum potential energy (QPE) 
in this metric. To date, in the literature, QPE has been attributed uniquely to microparticles. 
It has been confirmed that QPE would reduce to a time dilation in spacetime. For instance, 
when negative muons are captured in atomic s-states, their lifetimes are increased by a time 
dilation factor corresponding to the Bohr velocity. So, we expect that, if this QPE-induced 
time dilation is extended to macroscopic matter at cosmic scales, it may be disguised as 
some extra gravitational potential well. We will show that this extra time dilation can 
successfully explain the cosmic redshift and galactic rotation curve problem. For micropar- 
ticles, the troubles encountered all arise from the fact that when the QPE is significant, its 
gravitational effect is not geometric (i.e., mass-dependent). This results in the concept of 
the geodesic for test particles becoming ambiguous. Therefore, as mentioned earlier, if we 
want to apply this metric to gravitational systems composed of macroscopic particles, so 
that QPE is significant only at cosmic scales and can be neglected locally, we must establish 
a mechanism that connects the quantum effects of macroscopic matter to the spatial scales 
of gravitational systems. The mechanism, however, can be clearly demonstrated in the 
applications of our theory to cosmology and galaxies, as illustrated in the following section. 


3. Applications: Cosmology and Galaxies 


To investigate the mechanism that connects the quantum effects of macroscopic matter 
to the spatial scales of gravitational systems, we begin by considering a radially symmetric 
and static gravitational system consisting of N free-falling particles, each with an identical 
mass of m. We assume that the direct collisions and the close encounters due to gravity 
between the particles can be neglected. This allows us to approximate the gravitational 
effect of the particles by a smooth distribution of matter. Namely, we assume a smooth 
mass density function pm(r) and a smooth number density function e(r) defined by 


N 


Pm(r) = mp(r) = PG — Tri). (47) 


It is evident that the accuracy of the approximation improves as the mass m decreases 
and the particle number N increases. However, for a fixed mass density 0,,(r), reducing 
the mass m leads to an increase in the number density p(r). On the other hand, the scale- 
dependence of a particle’s mass in a gravitational system can be easily demonstrated. 
A particle with a mass of m contributes to the average mass density of a system with a size 
of À through the following relation: 


o(m,A) ~ ua (48) 


Clearly, from the perspective of average mass density, the “effective mass” of a particle 
decreases as the size of the system increases, following a trend of ~A~°. While this may 
seem trivial, it becomes crucial in our understanding of quantum effects at cosmic scales. 
For example, in our Milky Way, A~10 kpc. The Sun (m = Mo) contributes to the average 
density of the Milky Way as an electron (m = 9.1 x 10~78g) contributes to a system of size 
A~230 cm, a macroscopic size. So, as an assumption, we extend the quantum nature of 
electrons to macroscopic matter particles distributed in sufficiently large systems. 

In a gravitational system, the Dirac equation for a free-falling particle describes a 
particle of fixed proper mass m, a constant spin /2, and a proper probability density 
p(r). We consider the probability density as the number density of the system. In the 
case of spin-4 microparticles, the mass m can vary from neutrinos (m ~ 0) to neutrons, 
while the spin remains constant at h/2. Therefore, the spin of a particle is independent 
of its mass. This fact can be naturally extended to macroscopic particles. Typically, when 
dealing with large masses (macroscopic matter), the negligible quantum effect can be 
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understood as the short wavelength approximation, commonly referred to as the classical 
limit. In the context of radially symmetric gravitational systems, we can naturally interpret 
the negligible quantum effect as a result of the negligibly small number density of particles 
(since each particle has exactly a fixed value h/2 of spin). Given that reducing the mass m 
results in an increase in the number density p(r), we can introduce a decreasing function 
m(A) (A is the size of the gravitating system), which in turn leads to an increasing function 
p(r). Consequently, the strength of the quantum effect would increase as the distance (or 
scale) to an observer increases, aligning with our postulation. 

To be more precise, we define the density of QPE from (38) as (note that S? = —16?, 
and p(r) = pm(r)/m(A)) 


—2x( Ont)’ (49) 


where «(A) is a dimensionless parameter that depends on the scale A. We require that 
the strength of the QPE vanishes when A — 0 and increases with increasing A. Note 
that, in order to separate the size A dependence of pg completely from coordinate r, we 
express po(r, À) as proportional to p7,(r), and the constant mass density p; of the entire 
universe is introduced into the definition of og only for dimensional consistency. This 
choice has nothing to do with our contemplation of the notion that the entire universe 
may exert a cosmological influence on local galaxies, as elaborated upon in subsequent 
sections. Needless to say, the functional form of «(A) must be universal, meaning that 
it is independent of any particular gravitational system, although it can be different for 
the entire universe and galaxies. The specific properties of «(A) should be determined 
theoretically from first principles. However, as a first attempt to compare our predicted 
results with observations, when theoretical formulas for «(A) are lacking, we can derive a 
phenomenological formula by fitting the data to observations. 
From (38), we observe that the expression for M can be rewritten as 


M(x) =á | (om0) + pal’, A) rar 


r r 
=4r | Pm(r')r?dr' + se o2, (11? dr' (50) 
0 0 


c 


=Mnm (r) + Mo (r), 


where we have set A = r in the integrand to reflect the fact that, in PQ (r’, A), the scale A isr, 
the size of the “system” with mass M(r). And 


Malt) = an [ pnt") (51) 
4 g la ah f 
Mo(r) = sain [ Pim (1! )r'?dr (62) 


represent the conventional mass (energy) in general relativity and the QPE within r, re- 
spectively. From (49) and (50), it can be clearly seen that the distribution of QPE exhibits 
no extra symmetries, i.e., it aligns with the mass distribution. The reason for this is that 
the scale A characterizes only the amount of QPE for macroscopic matter at cosmic scales. 
Therefore, as the scale increases, the amount of QPE also increases, causing the quantum 
effects to vanish locally. This is in contrast to electrons, whose QPE can never vanish at any 
scale, thus resulting in possible anisotropy if not averaged [31,32]. 
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As previously mentioned, the functional form of «(A) should be determined based 
on observations. Nonetheless, we have found that the following general form is a suit- 
able choice: Al —1) 

e pen 
a (x ) = x2 
where Ho represents the Hubble constant, which will be defined subsequently. We leave the 
coefficients A, B, and C, possibly scale dependent, as free parameters to fit the observational 


data of the gravitational systems under consideration. 


+ B(e*—1)+C, with x = HA, (53) 


3.1. Cosmology 


Now we are ready to study cosmology. Assuming a static universe with a homoge- 
neous and isotropic matter distribution at large scales, characterized by a constant mass 
density pc, the mass M in (50) then becomes 


J r 
M(r) =4n | pcr’ dr’ +4ra(r) f per’?dr' 
0 0 


i (54) 
=5 ar(1+a(r)), 
where 8 
TT 
Hy = e: (55) 


is the Hubble constant. In cosmology, the entire universe can be regarded as an gravitational 
system with an arbitrarily large scale. To compare observations, we suggest setting the 


parameters in (53) as A = 2, B = 0, and C = —1. Namely, for cosmology, we assume 
2(eHor — 1 
ac(A) = a -1, (56) 


which yields desirable results: 
M(r) = r(e™" — 1) = Hor? + 5H Aes tg (57) 


Before formulating the metric for the spacetime of the universe, it is important to highlight 
that the term 4H47° resulting from pe in M(r) given in (54) should not be disregarded. One 
could argue that the universe is a distinct gravitational system, exhibiting a homogeneous 
and isotropic distribution of matter, and infinite in size. Consequently, the observation 
of gravitational redshift may not be possible under these circumstances. However, this 
conclusion is derived from Newtonian theory or general relativity, and does not hold true 
when quantum effects come into play. In fact, the term vetted can be included as the QPE if 
we choose C = 0 in (53). In any case, we should keep in mind that any choice of «-(A) must 
be subjected to testing against observations. Moreover, due to the presence of quantum 
effects at cosmic scales, we have to abandon the idea of the absolute spacetime predicted by 
general relativity. This is quite similar to the case of special relativity, when Einstein had to 
relinquish the concept of Newtonian absolute space and time. As a result, the conventional 
notion of the geometry of the entire universe becomes meaningless. For instance, it loses 
significance to classify the universe as open, flat, or closed according to its mass density pe. 
Thus, for cosmology, we have from (45) 


1/2 
= (1 7 ee) = @- 2eHor 1/2, (58) 


Hence, the metric (46) for the universe becomes 
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dt? =(3 — 2eo") dt? — (3 — 2e™0") 1dr? 


(59) 
— 1? (d8? + sin? 6d¢”). 

From this metric, it can be readily deduced that a horizon exists in the static universe, 
with the distance to any observer in the universe given by 


1 0.4 
n= In(3/2) ~ Ho’ (60) 
which is less than half a Hubble radius 1/ Hp. 

As a crucial outcome of our theory, we now derive the Hubble redshift for a static 
universe model. Suppose a light signal emitted from a source at r and received by an 
observer at r = 0, the redshift arising from the time dilation can be readily derived as 


vp £180) 1 
1+z A <= (61) 


It is evident that the redshift approaches infinity as the light source approaches the horizon. 
In the scenario where the source is located near an observer, the redshift can be approx- 
imated by a linear function of the distance r when z < 1, consistent with observations. 
To see this, we expand the expression in (61) as follows: 


z= Hor + 2H +.. (62) 


This shows that the first term, Hor, dominates the redshift at low values of z. In fact, it is 
precisely this favorable result that motivated our selection of the function «e (A) as assumed 
in (56). 

The close sources’ redshift z and distance r exhibit a well-established and fundamental 
observational fact of a linear relationship, one that remains independent of cosmological 
models. Hence, it is imperative to subject any plausible cosmological model to rigorous 
testing using relevant observations. We validate our cosmological model by comparing the 
predicted luminosity distance dy}, with the values derived from SN Ia data of the MLCS2k2 
Full Sample [36]. 

The luminosity distance is defined such that the Euclidean inverse-square law for 
the diminution in light with distance from a point source is preserved. Let L denote the 
absolute luminosity of a source at distance r and / denote the observed apparent luminosity, 
and the luminosity distance is defined as 


L 1/2 
E 9 


The area of a spherical surface centered on the source and passing through the observer is 
just 4mr2. The photons emitted by the source arrive at this surface having been redshifted 
by the quantum effect by a factor 1 + z. We therefore find 


E 1 
T In 12 = 
from which 
di, = r(1 +z)". (65) 


The coordinate distance r can be expressed in terms of z and Ho from (61), we thus obtain 


T "La- ap)] 


a= a, 2 1+2z)2 


(66) 
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The distance modulus is defined by 


p= mM =Siog( <S-) + 25, (67) 


where m and M are the apparent magnitude and absolute magnitude related to / and L, 
respectively, and we have explicitly specified that di, is measured in units of megaparsecs. 

We use the full cosmological sample of SNe Ia, MLCS2k2, presented in Table 5 of 
Riess et al.’s paper [36]. The sample consists of 11 SNe with redshift ranges from 0.01 to 0.1, 
which is ideal for testing our cosmological model. The only free parameter in our model 
is the Hubble constant Ho, with the best-fitted value being Hy = 59.0} 5 km/s/Mpc, 
as presented in Figure 1. 
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Figure 1. The fitting results of the predicted distance moduli given by (67) to the data of the MLCS2k2 
full sample [36]. The best fitted value of Ho is Hy = 59.0188, and x? = 0.0018. 


Our cosmological model fits remarkably well with the SNe Ia data, providing evidence 
that the cosmological redshift may be due to QPE-induced time dilation rather than the 
expansion of the universe. Our model is based on a static universe, which is made possible 
by the assumption of quantum effects at cosmic scales, as described by QPE, acting as a 
repulsive force to balance gravity. When Einstein introduced the cosmological constant A 
to balance gravity and obtain a static universe model, he failed, because A, regardless 
of its physical significance, is independent of the matter distribution. Consequently, any 
perturbations in matter result in unstable solutions of Einstein’s equations. In contrast, 
according to our proposal, QPE is not independent of matter distribution but is determined 
by the matter density, as shown in Equation (49). As a result, any perturbations in matter 
distribution result in corresponding perturbations in QPE, to balance the extra gravity, 
allowing the universe to remain static. 

It is worth noting from (49) that QPE depends on two independent factors: p%,(r) 
and «(A). As mentioned earlier, it is essential for «(A) to satisfy the condition that «(A) 
approaches zero as À approaches zero. Thus, even though galaxies possess a high but 
finite central matter density, the QPE would vanish in that region. However, in the con- 
text of cosmology, the QPE that is dependent on A displays some intricate yet justifiable 
characteristics. To see this, let us rephrase the function «,(A) in (56) as follows: 


1 
+H +- (68) 


2 
&e(A) = Hoà 3 
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It is evident that the first term 2/ HoA plays a crucial role in generating a linear correlation 
between redshift and distance at a low redshift. At first sight, the inverse relationship 
between this term and «,(A) may appear to contradict the requirement that «(A) approaches 
zero as A approaches zero. However, upon further examination, this does not turn out to be 
the case. To grasp this paradox, it is essential to recognize that it is Mg(A), not pg(r,A), that 
determines the metric and thus the geometry of spacetime for a specific observer, as shown 
in (59). 

In cosmology, the redshift of light signals emitted from distant sources is entirely 
attributed to the QPE, which varies with distance. However, this distance-dependent 
phenomenon does not imply any preferred direction in the universe, as demonstrated both 
theoretically and observationally. 

Although we prefer the presented static universe model, which lacks expansion, a big 
bang, and evolution of the entire universe, all local gravitational systems, such as galaxies 
and clusters of galaxies, can still evolve. In particular, for macroscopic systems, such as 
the solar system, planetary systems, and black holes, the evolution will continue to be 
governed by general relativity as usual. 


3.2. Our Solution to the Galactic Rotation Curve Problem 


After establishing the geometry of the entire universe’s spacetime, we now shift 
our focus to examining the galaxies. Traditionally, based on Newtonian theory, when 
a stellar system achieves a state of equilibrium, the kinetic energy supports the gravity, 
and this equilibrium state must adhere to the virial theorem. However, astronomical 
observations indicate that the gravitational force exerted by ordinary matter is insufficient 
to counterbalance the kinetic energy of the system. As a result, the presence of a mysterious 
substance known as dark matter has been postulated to account for the discrepancy in 
mass, often referred to as the missing mass problem. 

Based on the theory presented in this paper, we can propose an alternative explana- 
tion that challenges the need for dark matter. In fact, in the mass expression provided 
in Equation (38), the term QPE defined by Equation (49) serves as the counterpart to 
dark matter. 

The mass expression presented in Equation (50) is universally applicable to any radially 
symmetric system, making it suitable for investigating spherical galaxies. As mentioned 
earlier, «(A) should also be a universal function of scale A. We assume a functional form 
given by (53). In the context of cosmology, the scale can be regarded as arbitrarily large up to 
Th, and the reduced specific form of «-(A) provided in Equation (56) has been demonstrated 
to align with observational data. For galaxies, we use the parameters in (53) with A = 0 
and C = 0 to obtain 

&g(A) = B(e* — 1), (69) 


where, in this paper, we consider B as a constant, although it may be a parameter, possibly 
scale-dependent, to be determined by observations. When studying galaxies, the scale A of 
interest is much smaller than the Hubble radius 1/ Ho, thus «g (A) can be approximated by 
&g (A) = BHA. 

It is essential to account for the impact of the entire universe on local galaxies of 
finite mass, which is commonly referred to as the boundary condition when the metric of 
spacetime is concerned. Specifically, we mandate that the metric surrounding a local galaxy 
adheres to the condition that, as r approaches infinity, it converges to the metric of the 
entire universe (59), instead of the metric for flat spacetime, which has been conventionally 
used. To achieve this, we replace pm(r) in Equation (50) with pm(r) + pc. Since the impact 
of cosmological effects on local galaxies is unclear, we make the assumption that 


A(A)(eHo4 — 1) 
Hp? 


a (A) = (70) 
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when considering pc, allowing us to determine the function A(A) based on observations of 
galaxies. We thus have from (50) 


r r Amtae(r) øT 
M(r) =4n f intr’) + pe)rdr’ + 4reac(r) | per dr + ne ) f o2 (rr ?dr' 


: (71) 
=M,(r) + Mn (r) + Mo(r), 
where 
a ee 
Mn (r) =47 f Om (r rdr’, 
or _ 72 
vay BD Pay yar = 


(eHor — 
a 1) df ee) sa! 


Evidently, M,(r) signifies the contribution of the quantum effects of the background mass 
density of the universe, Mn (r) denotes the conventional mass, and Mo (r) represents the 
quantum effects of a galaxy itself. We thus have 


“mE o o 
(73) 


2Mn(r) Mol) = 


=[1— ater —1) 


We require gı for galaxies to satisfy the boundary condition, i.e., when r — œ, 
g1 — (3 —2e%0")1/2 the value for cosmology as shown in (58). For galaxies of finite mass or 
finite size, the last two terms ein) Mat > 0 when r — œ. We thus require A (r) — 2 
when r —> œ. This condition imposes a nature constraint on the cosmological effect on local 
galaxies, as shown subsequently. 

Our findings indicate that, from the perspective of mass-energy, the contributions of 
QPE to the metric exhibit no observable distinctions from conventional mass. Both QPE 
and conventional mass contribute to time dilation and distance contraction in precisely an 
identical manner, as demonstrated in Equation (46). As such, it is reasonable to employ the 
conventional approach from general relativity when considering the stress-energy tensor 
and geometry of spacetime. 

Although our findings are presented in relativistic form, the transition to a non- 
relativistic scenario for galaxies is straightforward. Actually, we can directly utilize the mass 
expression provided in Equation (71) and apply Newtonian theory to investigate galaxies. 

Obviously, near the center of a galaxy, the gravitational effects arising from Mp (r) and 
Mo (r) can be neglected when compared to those of Mm(r). Towards the outer regions, 
both M;,(r) and Mo(r) become significant, while the contribution of M,,(r) diminishes 
with increasing radius r. 

It is interesting to investigate the circular velocity of a test particle around a galaxy, 
which can be approximated by 


M(r) 
A(r) 
2 


v(r) = 
B327? fT > (rar! 79 


4 r 
Hor 4 A Pm(r)r?dr' + 
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where we have neglected the terms with Hjr" when n > 2 for the scale of galaxies. For a 
galaxy with a finite size a, when r increases to the outer regions with r >> a, we have 


v(r) = T + Mg(a)/a. (75) 


This expression does not imply a constant value of v(r) for large radii, as the Hor term 
in (75) becomes more and more significant with increasing r. In particular, this leads to a 
universal centrifugal acceleration (recall that A(r) — 2 for r — oo) 


2 
Z = cHo, (76) 


where c is the speed of light. This result reflects the fact that the matter in the entire universe 
can have a local observable effect in galaxies, a fact first discovered in MOND theory but 
that cannot be explained within the theory itself [10,37]. In fact, our cosmological model 
yields the value of Hp = 59.0, resulting in cHọ = 5.74 x 1078cm/s?. A value which is 
close to, but obviously larger, than that found for the universal acceleration parameter 
ay = 1.2 x 1078 cm/s? of the MOND theory [9,38]. This fact has a natural explanation. 
As an increasing function of r, A(r) — 2 only at cosmological distances, typically hundreds 
of Mpc. The value of ap is smaller than cHo only because it is detected at distances in the 
outer regions of galaxies, which are much smaller than the cosmological distance. Notably, 
conformal gravity yields comparable results when fitting rotation curve data [13,39]. In our 
theory, these phenomena are not considered mysteries. The cosmological effects on local 
galaxies are the natural results of the boundary condition, which is determined by the 
requirement that the spacetime metric around an isolated galaxy, created by the mass 
(energy) specified in (72), should align with that of the entire universe, as demonstrated 
in (59), when r — œ. 

Despite these remarkable achievements, the most effective way to validate our theory is 
by leveraging the abundant observational data on the rotation curves of flattened dwarf and 
spiral galaxies. As our results were obtained for spherical systems, they must be converted 
into axisymmetric systems where cylindrical coordinates are more suitable. One could 
replicate the procedure for axisymmetric systems [21,40], similarly to what was performed 
for spherical systems in this paper. However, a more efficient approach to achieve our 
objectives is to consider the relevant terms on the right-hand side of Equation (74) as the 
gravitational potential generated by a point mass (or a mass element). Subsequently, we 
calculate the total potential for disk galaxies in cylindrical systems (R, Q, z) [13]. 

The first term on the right-hand side of (74) is A(r)cHor/2. This is a linear potential 
originating from the quantum effect of the entire universe, and thus is independent of 
specific local galaxies. We simply express the contribution of this term to the total circular 
velocity of a test particle on a thin-disk as 


v? (R) = A(R)cHoR. (77) 


In this paper, we initially overlook the gradually increasing nature of the function A(R) 
and consider it as a constant Ag, which can be determined using rotation curve data. 

The second term on the right-hand side of (74) is the usual Newtonian potential. We 
write the corresponding potential in cylindrical coordinates as 


co 27 oo R’o (R’ z’) 
Dn (R, =-G f ar’ f ag f dz! mAy , (78 
m(R,2) 0 ae Oe (R2 + R’? — 2RR' cos p! + (z — 2!)2)1/2 oe 


where G is the gravitational constant (we temporarily transition back from natural units 
from now on in this subsection). Inserting the cylindrical coordinate Green’s function 
Bessel function expansion 
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1 
(R2 + R’? — 2RR' cos! + (z — 2')2)1/2 


=) [ dk Jn (KR) Jn (KR) #9) -k=l (79) 


into (78) yields 


Dy (R,z) = —2"G f dk J dR’ I dz! R'om(R', 2) Jo(kR)Jo(kR')e =l. (80) 
0 0 —oo 


For razor-thin exponential disks with pm(R,z) = =(R)d(z) = Xpe~®/®05(z), where Xo is 
the central surface mass density and Ro is the disk scale length, the potential is 


Pin(R) = —7GLoRoy[Io(y/2) Ki (y/2) — h(y/2)Ko(y/2)], (81) 


where y = & , In and K, are modified Bessel functions. If we differentiate Equation (81) 
with respect to R, we obtain the circular velocity contributed by the Newtonian potential 


2 
chy) = REM = FE (SES) (FE R Uto(u/2)Koly/2) — (y/2)Ki(y/2)]- 2) 


Another potential of the quantum effect for galaxies on the right-hand side of (74) is 
(6477/3Hp) Jj Pm (r')r?dr'. This corresponds to a logarithmic potential, and the circular 
velocity for the razor-thin exponential disk is 


B87G? 
ral 2(R 
Hake ae [a rf” ap iz aeRO) 


(83) 
_ B4r*G?X3Ro 1- (14.28) 28/20], 
3cHo Ro 
The total circular velocity used to fit the data is therefore 
o(R) = \/02(R) + 03,(R) + 03,(R). (84) 


We fit the predicted circular velocity to the data provided in the Spitzer Photometry 
and Accurate Rotation Curves (SPARC) database [41]. The SPARC database is the largest 
sample, to date, of rotation curves for every galaxy. It isa sample of 175 nearby galaxies with 
new surface photometry at 3.6 mu, and high-quality rotation curves from previous HI/H 
studies. As a first try, we consider each of the 175 galaxies in the sample as a razor-thin 
exponential disk characterized by a central surface mass density Xp and a disk scale length 
Ro. In addition to Xp and Ro, we treat A(r) = Ag in (77) and B in (83) as free parameters in 
our fitting. Figure 2 displays 12 of them, including 6 low-surface-brightness (LSB) galaxies 
(6 upper panels) and 6 high-surface-brightness (HSB) galaxies (6 lower panels). 

The fittings reveal that Ag = 0.02 for all galaxies, with B = 3.0 for LSB galaxies 
and B = 0.3 for HSB galaxies. The cosmological effect on local galaxies thus suggests 
an acceleration of AgcHy~10~° cm/s?, significantly smaller than the universal value of 
cHy = 5.74 x 1078 cm/s?. Our fittings further suggest that Ap can be substituted with 


Ao + 2(R/Mpc)? 
1+ (R/Mpc)2 


A(R) = (85) 
This slowly varying function provides a good fit to the rotation curve data and also ful- 
fills the requirement that A(R) converges to 2 as R approaches cosmological distances. 
On the other hand, the value of B for LSB galaxies exceeds that for HSB galaxies, sug- 
gesting that LSB galaxies require relatively more QPE than HSB galaxies, and B(A) is also 
scale dependent. 

All these results indicate that, in terms of mass-energy, scale-dependent quantum 
phenomena do not exhibit a linear superposition, but rather a hierarchical process. 
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Figure 2. The best-fitting rotation curves of LSB (upper 6 panels) and HSB (lower 6 panels) galaxies 
with the model presented in this paper. The panels are listed by the increasing effective surface 
brightness of the galaxies. In each panel, the dotted curve shows the contribution from the cosmolog- 
ical quantum effect, given by (77); the dashed curve indicates the contribution from the luminous 
Newtonian potential given by (82); the dash-dotted curve shows the contribution from the quantum 
effect of the galaxies themselves given by (83); and finally, the solid curve is the total circular velocity 
given by (84). 


Despite the very simple model used to describe the matter distribution in LSB and 
HSB galaxies, the fittings can still capture the main phenomenological results observed in 
other theories of gravity. For instance, the quantum effect contributions (v? + Vo) dominate 
most regions of LSB galaxies, whereas this is only true for the outer regions of HSB galaxies. 

Our proposed quantum effect on large scales has proven to be successful in explaining 
the mass discrepancy problem in galaxies. 


4. Conclusions and Discussions 


We have proposed the quantum nature of spacetime at cosmic scales. To investi- 
gate this postulation, we examined spherically symmetric and static gravitational systems 
composed of free-falling particles with identical mass m. We assumed that the matter dis- 
tribution can be smoothed using the proper number density o(r) and proper mass density 
Pm(r). We further assumed that the Dirac spinor y(r) can fully capture all the physical 
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aspects of the particles and satisfies the covariant Dirac equation. Therefore, we recog- 
nize p(r) = P(r) P(r) and pm(r) = mp(r)p(r) = mp(r). For a given mass density pm(r), 
the interplay between the number density p(r) and the particle mass m offers us a inherent 
mechanism to establish a connection between the spin density S(r) = Lu(r)iysp(r) and 
the scale A of gravitational systems. 

We utilized gauge theory gravity (GTG) [26,27] to achieve the subsequent derivations. 
After solving Einstein equations, the metric of the spacetime for the radially symmetric 
and static gravitational systems was derived, as shown in (46). The fundamental results 
obtained thus far require additional physical interpretations and practical applications. We 
found that the quantum potential energy (QPE) defined by (49) 


polr,A) = ~5xS%(r) = ph(r)a(A)/pe 


is contained in the metric via the mass (energy) M(r), where «(A) is a dimensionless 
function of the scale of the system involved. The mass within radius r can then be expressed 
in (50). 

It is important to note that at each radius r, we consider M(r) as the mass of a 
gravitational system; therefore, we identify r as the scale A of the “system”. As such, 
for any given mass density pm(r) of a gravitational system, we can derive the metric that 
encompasses not only the classical mass but also the quantum potential energy (QPE). 
When the scale of the gravitational systems are macroscopic, for instance, the solar system, 
the QPE can be neglected, and the Einstein-Newtonian theory of gravity is recovered. 

The dimensionless function «(A) is universal, and it should be determined theoretically 
from first principles. However, when theoretical formulas are lacking, we can derive a 
phenomenological formula by fitting data to observations. We discovered that the general 
expression presented in (53) successfully achieved our objective. 

While our fundamental results were initially derived for radially symmetric and static 
systems, they can be applied to any static gravitational system. By substituting the proper 
mass density pm with Om + po + +c + Pcltc(A), where pe is the constant mass density of 
the universe, we establish a fundamental boundary condition that the spacetime metric 
encompassing any local gravitational system must adhere to. Remarkably, this boundary 
condition inherently gives rise to the cosmological impact on local galaxies, as evidenced 
by observations of galactic rotation curves [9,10,13,38,41,42]. 

When applied to cosmology, our model yields a static universe. The predicted luminos- 
ity distance—redshift relation fit remarkably well with SNe Ia data, providing evidence that 
the cosmological redshift may be due to QPE-induced time dilation rather than expansion 
of the universe. It is not surprising that when we extend Dirac’s theory for free electrons to 
macroscopic particles (including stars within galaxies and galaxies within the entire uni- 
verse), we are effectively moving away from the notion of a preferred absolute spacetime as 
suggested by general relativity and commonly accepted by many physicists, as illustrated 
in the standard ACDM cosmology. Without the existence of absolute spacetime, physical 
processes should be described from the perspective of any observer in the universe inde- 
pendently and equivalently, without relying on the absolute spacetime reference frame of 
the universe. Therefore, the state of a celestial body should be characterized in the frame 
associated with a specific observer rather than in the frame of the absolute spacetime of 
the universe. Consequently, we assume that, for any observer in the universe, a distant 
celestial body behaves akin to a micro-particle and can be effectively described using 
Dirac’s theory. On the contrary, in a scenario where absolute spacetime exists, all matter 
particles, including observers, must be described in relation to it. As a result, the motion 
of distant celestial bodies is typically explained by a combination of the local velocity and 
the overall expansion velocity. In this context, we should not anticipate quantum effects at 
large scales for any observers. 

When considering galaxies, it is intriguing to note that quantum effects can serve as 
a substitute for dark matter. To validate our theory, the most effective approach was to 
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utilize the extensive observational data on the rotation curves of dwarf and spiral galaxies. 
Remarkably, our proposed quantum effect at large scales successfully addressed the mass 
discrepancy issue in galaxies. Specifically, our theory presented a fitting formula, as 
outlined in (85), which elegantly explains the correlation between the universal acceleration 
cHy suggested by cosmology, 49 proposed by MOND, and 0.02cHp resulting from the 
cosmological effect on local galaxies [13]. 

Although not investigated in the present study, it is interesting to note that there is 
another very important piece of evidence supporting our new theory about the quantum 
effects of macroscopic matter at cosmic scales. This evidence comes from a large-scale 
structural survey of the universe, which reveals the fractal geometry of matter distribu- 
tion at these scales, as cited in Ref. [43]. It has now been established that the quantum 
randomness of the electron can be mimicked by Brownian motion; the simulated trajectory 
is continuous but non-differentiable. In fact, the non-differentiable trajectory or path of an 
electron exhibits a fractal structure due to the uncertainty principle between its position and 
momentum in the conventional matrix and tensor versions of quantum mechanics [44]. Nat- 
urally, the observations of the fractal geometry of the large-scale structure of the universe 
can be regarded as independent and compelling evidence that supports our assumption 
when extending the reasoning about the relationship between the quantum nature of matter 
and its fractal path from microscopic matter to macroscopic matter. 

In conclusion, we have expanded the quantum characteristics of microparticles to 
encompass macroscopic matter at cosmic scales. As illustrated, this expansion was remark- 
ably coherent. The subsequent outcomes were derived within the framework of the firmly 
established theory of gravity (GTG) and mathematical language (STA). While there is room 
for improvement in the technical details of our derivations, the fundamental discoveries 
are credible and, notably, were verified through astronomical observations of SNe Ia and 
galactic rotation curve data. In particular, our proposal provides an alternative view point 
to understand the quantum nature of spacetime. 

The manifestation of physical laws should be independent of the mathematical lan- 
guage used, and thus we believe that our proposal could also be accomplished, for instance, 
through the tensor analysis approach. Of course, using an unfamiliar language may make 
it challenging for readers to comprehend our theory. However, it has been shown that STA 
and GC can significantly streamline the calculations in our area of concern. Additionally, 
STA can illustrate the cohesive link between classical mechanics and quantum mechan- 
ics, thereby enhancing our comprehension of the enigmatic aspects of quantum theory 
and spacetime. 

For future work, we could test our theory through observations of gravitational 
lensing, cosmic microwave background radiation, and more sophisticated models for 
rotation curves. In particular, we could investigate the bullet cluster problem, which is 
considered the best evidence for the existence of dark matter. A more ambitious project 
would be to derive the «(A) formula from first principles. 
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Appendix A. Spacetime Algebra and Geometric Calculus 


Spacetime algebra (STA) is generated by a 4-dimensional Minkowski vector space. 
The inner and outer products of the four orthonormal basis vectors {Y,,/ = 0...3} are 
defined to be 


1 ; 
Yu: v Z 5 uY + Yn) = Nw = diag(+ — —-) 
1 
WAW = z MY = YY): (A1) 
A full basis for the STA is 
1, {Yu} {kik}, {ifu} i (A2) 


where op = YkYyo,k = 1...3 constitute the orthonormal basis of a 3-dimensional vector 
space relative to yo, and i = 70717273 = 710203 is called pseudoscalar. The STA is a linear 
space of dimensions 16. We refer to the general elements of STA as multivectors, and each 
multivector decomposes into a sum of elements of different grades. If a multivector contains 
only grade-r components, we call it homogeneous, and is denoted by A;. We call grade-0 
multivectors scalars, grade-1 vectors, grade-2 bivectors and grade-3 trivectors. A grade-r 
multivector is called simple or an r-blade if and only if it can be written as 


1 
Ay = ay A ag A... Alar = JI Li(senm)an(ay@n(a)-An(r) (A3) 
where 7 is a permutation of 1 through r, ‘sgn 7’ is the sign of the permutation, 1 for even 
and —1 for odd, and the sum is over all r! possible permutations. For r = 2, this reduces 
the outer product of two vectors 


a/b = ~(ab—ba), (A4) 


NIe 


which is in agreement with the definition in (A1), and clearly a ^ a = 0. 
The geometric product of a grade-r multivector A; with a grade-s multivector Bs is 
defined simply by A;B;, which decomposes into 


A,B; = (A;Bs) rs] + (A;Bs)\p—s|42 Hart (ArBs)r+s, (A5) 


where (X), denotes the projection onto the grade-r part of X. The grade-0 (scalar) part 

of X is written as (X). We employ “:” and “^A” symbols to denote the lowest-grade and 
highest-grade terms in (A5), so that 

Ar $ Bs = (ArBs) |r—s|7 (A6) 

Ar ^B; = (ArBs)r+s, (A7) 


which are called inner and outer products, respectively. They represent a generalization of 
the geometric product between two vectors a and b, which can be expressed as follows: 


ab=a.-b+a^b. (A8) 


We define the reverse of a geometric product AB by (AB)~ = BA, so that for vectors 
41,42, ..., Ay, we have 
(a1 42...ay)~ = Arap—1..-41. (A9) 


This reverse operator obeys the rule 


(A+B)~ = A+B, (A10) 
(A) = (A), (A11) 
i = a, (A12) 
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where a is a vector. It is easy to show that 
A= (18s. (A13) 
Thus, suppose r < s, the inner and outer products satisfy the symmetry properties 


A,:Bs = (—-1)'6-)B,-A,, (A14) 
A; A Bs (1° BA Ar. (A15) 


The scalar product is defined by 
A x B= (AB). (A16) 


From (A6), A, * Bs is nonzero only if r = s, thus the scalar product (A16) is commutative 
(AB) = (BA), (A17) 


which proves to be very useful in our calculations. The inner product and outer product 
obey the distributive rule 

A- (B+C)= = A-B+A-C (A18) 

AA(B+C) = AAB+AAC. (A19) 


The outer product is associative 
AA (BAC) =(AAB)AC. (A20) 
The inner product is not associative, but homogeneous multivectors obey 


Ay-(Bs-C:) = (Ar^ Bs): Ci for r+s<t, (A21) 
Ar: (Bs-C:) = (Ar: Bs): Ci, for r+t<s. (A22) 


For any vector a, it can be proved that 


1 
a- Ar = GA) 1 = 5 Ar — (—1)"A,a), (A23) 
1 
aN Ay = (aAr)r1 = 5 Ar + (—1)'A;a). (A24) 
From which, we immediately have 
aA, =a. A, +a ^A. (A25) 
We also have the following useful identities 
a: (A,B) = (a-A,)B+ (-1)'A;(a-B) (A26a) 
= (a^Ar)B-—(—1)"A,(a ^B) (A26b) 
A (A,B) = (a^Ar)B-—(—1)"A,(a- B) (A26c) 
= (a-A,)B+(—1)"A,(aA B). (A26d) 


These identities imply the following particularly useful relations: 


a- (Ai) = (aAA,)i, (A27) 
aN\(A;i) = (a-A,)i, (A28) 


where i is a pseudoscalar, and we have used the fact that a ^ i = 0. 
Any vector a can be decomposed in terms of {Yy} into 


a =A: Yuy" =a- Y yy, (A29) 
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where the summation convention is implied. Similarly, an arbitrary multivector A can be 
decomposed into 
A= L Apv YÉN NY, (A30) 
u<...<v 
where 
Apv =A: (Yv N A Yy). (A31) 


We further define the commutator product 


A x B = (AB — BA), (A32) 


NI =e 


which satisfies the Jacobi identity 
Ax(BxC)+Bx(CxA)+Cx(AxB)=0. (A33) 
In addition, we have the Leibnitz rule for the commutator product 
A x (BC) = (A x B)C+ B(A x C). (A34) 


This rule is particular useful with bivectors. It can be proved that 


A? x Ar = (A24; = A,A?) = (A2Ar)r, (A35) 


N| rR 


which means that if one of the factors is a bivector, then the commutator product preserves 
its grade. Thus, for (A); = 0, we have 


AzA =A- A+A X A+A2/AA. (A36) 


Since the commutator product with a bivector is grade preserving, the identity in (A34) still 
holds if A = Az and we replace all geometric products with either inner or outer products: 


Az X (BAC) = (Az x B)\C+BA (Az xC). 
Finally, for any vector a and trivector T, it can be proved that 
(a-T)xT=0. (A38) 


Geometric calculus (GC) is the extension of a geometric algebra (like STA) to include 
differentiation and integration. Let the multivector F be an arbitrary function of a multi- 
vector argument X, then the derivative of F (X) with respective to X in the A direction is 
defined by 

F(X + TA) — F(X) 


Ax*dxF(X) = lim = j (A39) 


where the multivector partial derivative 0x inherits the multivector properties of its argu- 
ment X. Suppose that the set {e,} form a vector frame (which is not necessarily orthonor- 
mal), the reciprocal frame is determined by [45] 


ef = (-1) te Neg A... AG Aw Ane! (A40) 


e =e] ^e N... Nen (A41) 


and the check on é denotes that this term is missing from the expression. As usual, the 
two frames are related by 
ej = OF. (A42) 
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From this frame, the multivector derivative dy in Equation (A39) is defined by 
dx= Vo eA... AeA... Ae) * Ox. (A43) 
1<...<j 
The Leibnitz rule can be written in the form 
dx(AB) = 0xAB + 0x AB, (A44) 
where the overdot indicates the scope of the multivector derivative. 
We have from (A39) and (A43) 
dx (XA) = Px(A), dx(XA) = Px(A), (A45) 


where Px(A) is the projection of A onto the grades contained in X. These results are 
combined using Leibnitz’s rule to give 


dx (XX) = dx (KK) +x (XX) = 2X. (A46) 
For a vector argument x and a constant vector a, (A39) and (A45) yield 
a-OxxX =a =0x(xX-a). (A47) 


For a vector variable a = aly, = a- Yuy" = a- "Vy, where y” constitutes the 
reciprocal basis and satisfies Yy -y" = Ons the vector derivative can be defined as 


d 


a= dal’ 


(A48) 


For the derivative with respect to a spacetime position vector x, we use the symbol 
V =ð = y gr, if x = xy. From (A47) and (A48), we can obtain useful results 


da = dpb Da = Y Yp da: (A49) 


Some results for the derivative with respect to position vector x in an n-dimensional space 
are [45] 


dx(x- Ar) = TA, (A50a) 
d(x A Ay) = (n-=-r)A, (A50b) 
xArt = (—-1)'(n—2r)Ay. (A50c) 


From (A21) and (A50a), an r-blade A; can be expressed as 
1 
Ar = heal A... A 0a, (ar A. A a1) * Ar. (A51) 


When considering a vector argument x and a constant vector a, (A39) becomes the 
definition of the directional derivative a- V, thus 


eVP=eeO< tee) Hie CE, 


A52 
dt t=0 170 T 02) 


Then, the general vector derivative can be obtained from the directional derivative us- 
ing (A49) as 

VF = 0xF(x) = Oga- OxF (x). (A53) 

The directional derivative (A52) produces from F a tensor field termed differential of 


F, denoted variously by 
f(a) =F, =a. VF. (A54) 
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The underbar notation serves to indicate that f(a) is a linear function of a. This induced lin- 
ear function is very important for us to describe the apparatus of GC for handling transfor- 
mations of spacetime and the induced transformations of multivector fields on spacetime. 
Suppose there is a diffeomorphism that transforms each point x in some region of 

spacetime into another point x’ as 
x = f(x). (A55) 


This induces a linear transformation of tangent vectors at x to tangent vectors at x’ given 
by the differential 
a'(x') = f(a) =a- Vf. (A56) 


If we regard x as a map representing the ordering of points in spacetime, then x’ can be 
interpreted as a different map, or a remapping of the same spacetime. The transformation 
f also induces an adjoint transformation f, which takes a tangent vector b’ at x’ back toa 
tangent vector b at x, as defined by 


b(x) = f(b’) = dxf (x) -b'(x’). (A57) 
The differential and its adjoint are related by 
b' -f(a) =a f(b’). (A58) 
By using this relation, we can find one transformation from another by 
f(a) = db- f(a) =d (f(b): a), (A59) 
f(a) = db- f(a) = (F(b) - a). (A60) 


In addition to the induced linear transformations f(a) and f(a) of tangent vectors, 
by the rule of direct substitution, (A55) can also induce a transformation of a multivector 
field F(x) defined by 

F'(x') = F'(f(x)) = F(x), (A61) 


in which directional derivatives of the two functions are related by the chain rule 


a:-VF = a-dxF'(f(x)) (A62) 
= (a: Vxf (x): Ov F(x’) (A63) 
= f(a)-V'F =a-f(V')F (A64) 
a’ VIFF. (A65) 
The operator identity is 
a-V = f(a): V =a-f(V')=a'-V (A66) 


Differentiation with respect to the vector a yields 
V = f(V') or V' =f" (¥). (A67) 


Now is an opportune moment to discuss linear algebra. In fact, GC enables us to 
carry out coordinate-free calculations in linear algebra, eliminating the need for matrices. 
Every linear transformation f on spacetime has a unique extension to a linear function on 
the whole STA, called outermorphism. For arbitrary multivectors A, B, and any scalar a, 
the outermorphism is defined by the property 


f(AAB) = f(A) A f(B), f(a) =a. (A68) 
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It follows that for an r-blade A; = a1 A ... A ay, 
f(Ar) = fai) A... A f (ar). (A69) 
Since the outermorphism preserves the outer product, it also preserves grade: 
f((A)r) = (F(A))r (A70) 
for any multivector A. This implies that f alters the pseudoscalar i only by a scalar multiple: 
f(i) = (detf)i, (A71) 


which defines the determinant of f. The product of two linear transformations h = gf also 
applies to their outermorphisms. It follows from (A71) that n 


h(i) = g(f(i)) = (detf)g(i) = (det f) (det g)i (A72) 
det(gf) = detgdetf (A73) 
detf™! = (detf)™'. (A74) 


Every f has an adjoint f, which can be extended to an outmorphism denoted by the 
same symbol 


(Af(B)) = (Bf(A)) (A75) 


for any multivectors A and B. Unlike the outer product, the inner product is not generally 
preserved by outermorphisms. However, it obeys the law 


Ar-f(Bs) = F(f(Ar) Bs) for r< s, (A76) 
f(Ar)-Bs = f(Ar-f(Bs)) for r >s. (A77) 


From these identities, we can construct a formula for the inverse of a linear transformation. 
Consider a multivector B, lying entirely in the algebra defined by the pseudoscalar i, 
we have 


f(i)B = (det f)iB = f (if (B)). (A78) 
Replacing iB by A, we find 
(det f)A = f(if(i-'A)) (A79) 
with a similar result holding for the adjoint. It follows immediately that 
f(A) = if@tA)(det f), (A80) 
f(A) = if(i*A)(detf)7?. (A81) 


Once again, one great advantage of GC is that it eliminates unnecessary conceptual 
barriers between classical, quantum, and relativistic physics. 


Appendix B. Stress—-Energy Tensor Derived from Dirac Theory 


In this Appendix, we show that the stress—energy tensor derived from the Dirac 
equation for free spin-5 particles can be decomposed into the sum of a symmetric part and 
an antisymmetric part. The symmetric part represents a classical pressureless ideal fluid, 
while the antisymmetric part represents the contribution of quantum potential energy. As a 
result, when spin is zero, the stress-energy tensor becomes that of a pressureless ideal 
fluid. We follow Hestenes’s work [18], except that we set = 0 for free spin-} particles. 
As mentioned in our previous paper, the spinor field at spacetime point x takes the form [14] 


p(x) = p(x)'/?R(x), (A82) 
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where p(x) is a scalar, representing the proper probability density and R(x) is a rotor 
(Lorentz rotation) satisfying RR = 1. The rotor R can be used to transform a fixed frame 
{Yu} into a new frame {eu} 

en = RY R. (A83) 
We identify v = eg as the proper velocity associated with the expected history x(t) of a 
particle, v = a Our present objective is to derive a stress—energy tensor for a pressureless 
ideal fluid from the Dirac theory. The desired form of the tensor is as follows: 


Tyv = Pm0pPv, (A84) 


where v = vy” and pm = mp, as explained in the previous paper. 
For later use, we define 


1 | 
spin bivector:S = aRrnk = 56201 (A85) 
1 g 1l 
spin density trivector : S = z PTY = 5 pe2e19 = pSv. (A86) 
The Dirac equation is 
Vpiy3s = my. (A87) 


A stress—energy tensor T(a) is a linear vector function of a vector variable a, which denotes 
a flux of energy-momentum through a hypersurface with normal a at spacetime point x. 
The tensor T(a) for the free Dirac field is given by [18,21] 


T(a) = y” (adv piysp) = ya- (vpiysh)1- (A88) 
This tensor is not symmetric, and its adjoint (transposed tensor) is 
T(a) = d(T (b)a) = (a- Vypiysh)1. (A89) 


We can decompose T (a) into a sum of a symmetric part and an anti-symmetric part as [45] 


T(a) = Ts(a) + Ta(a), (A90) 
where i i 

Ts(a) = 5(T(a) + F(a) = 58a(a-T(a)), (A91) 

Tala) = 5(T(a) — T(a)) = Fa- (r ATC)) = -2a (s ATO). (A92) 


From (A86) and (A89), we write 
da A Pla) = Y A (Hirsh 
= SV pirsh — 1" (pirs(Bnp)~))2 
= HV pis +7" (Yirs(@up)”))2 (A98) 


1 ae 
= —5(V(piy3))2 
=—V-S3, 
where we have used the fact that the Dirac Equation (A87) implies (Vpiy3t))2 = 0. Thus, 
according to (A92) and (A93), the anti-symmetric part of T(a) can be written as 


Ta(a) = 52-(V-S3) = 5(aAV) Ss. (A94) 
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We are now ready to define the local energy-momentum p, one of the most fundamen- 
tal quantities of the Dirac theory, as [18] 


pp =T(v) = oF T,, (A95) 


where v = vy, and T, = T(z) are understood. Now T, can be decomposed in the 
form [18] 
Ty = PVup + Nu, (A96) 


where N, = N (yp) describes the flow of energy momentum normal to the velocity. This 
can be verified from (A95) and (A96), which reads N(v) = N(v!7,) = v0! N, = 0, since 
voy = v2 =1. 

According to the conservation of energy and momentum, the divergence of a stress- 
energy tensor must be equal to the force density exerted on the particle. This relationship 
holds true even for free particles, where no external forces are acting on them. We thus write 


Vist Viet yo) So" = 0, (A97) 


where we have used the fact that the divergence of the anti-symmetric part of the stress— 
energy tensor vanishes. This is because, according to (A94), T, (V) =0. 

One advantage of STA is that the formulations of physical laws in Dirac theory, such as 
conservation laws and dynamics, are expressed in the same form as in classical mechanics. 
From (A83), it follows that 

a- Vey = O(a) X ey, (A98) 


where 
Q(a) = 2a-VRR (A99) 


is a bivector valued function of a vector variable, which can be explained as the angular 
velocity of the frame {e, } rotates in the a direction, and A x B = (AB — BA). From (A88) 
and (A96), we write 


Tuv = Ty: Ww = Yu: (dvpiys))1 = PPnPy + Nw, (A100) 


where py = p- Yu and Nyy = Ny - yv. To express T,,, with observables, we need to calculate 
d, piy3) from (2). We have 


dvpiysp = dv (P Riy R 


1 ow 5 
= zvPRiyaR + pd) Riy3R 


= (0 inp)S5 405) (A101) 
= (ð, Inp)pSv + pA, Sv 
= 0(W, -dS +0,S)v, 
Here, Oy = Q(yv) = 23 RR; W, represents a bivector defined as 
1 
W, = pee) (A102) 


which is called quantum potential, and this holds a pivotal position in quantum mechanics 
in the approach of causal interpretation [18]. The vector part of (A101) gives 


(dvpiysh)1 = Ov: S3 
= p(W,-v—%S-0+QO,-Sv+ (Oy x S)-v) (A103) 
=0(W,-v+Q,- Sv), 
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where we have used 0)S = Q, x S, which can be proved as follows 
1 5 
duS = z?n (RyrR) 
1 5 4 
= =((0,R R+R ðu R 
7 (( u Jy yik + RYN u ) (A104) 
1 
= On x S. 
We thus obtain from (A88), (A95) and (A103), 
T(v) =v- (OvWiysp)1y" = (v AQ»): S37” = pQ : SY” = pp, (A105) 
which means 
Pv = Oy -S. (A106) 
Finally, from (A100), (A103) and (A106), we obtain [18] 
Tyv = Vy py +P(A Yy): Wv. (A107) 


To achieve the desired result (A84), we need to identify the constraints that would align 
the momentum and velocity (i.e., v A p = 0), and ensure the symmetry of the stress-energy 
tensor (i.e., T (a) = T (a) or Tw = Typ). 

Let us first multiply the Dirac Equation (A87) on the right by ¢ to obtain 


Vy = mov. (A108) 


Next, we evaluate V yy2y1Ọ by substituting y from (2), resulting in 


 (Onbr2nh) = y (uV R210 + Ve(0nR) 121) 
= 7" ((Oup)S + pOA,S) (A109) 
= (Vp)S + py", S. 


The vector part of this equation gives 


(Vibyon)1 = ((Vpe)S + py OnS)1 


Vp): S+py#Oy,-S + py! (Qu xS 
(Ve) S +py"Qy: S+ py": (Oy x S) (A110) 
(Vp) -S+ppt+pv-S 
= V - (pS) + pp. 
From (A108) and (A110), we obtain 
p = mv — v - (pS). (A111) 


Decomposing the stress—energy tensor into its symmetric and anti-symmetric parts 
in terms of observables proves to be highly advantageous. In the main text, we specifi- 
cally utilize the adjoint form outlined in this Appendix, hence we express T,,, as derived 
from (A107) as follows: 

Ty = ppyv + [dp (S30)] - v. (A112) 


By substituting p given by (A111), we have 
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T(a) =pa- pv + [a - V(S3v)] -v 
=p|ma -v — (a A V): (S3v)|v + [a- V(S30)] -v (A113) 
=Pma-vv + [a - V (S3v)] -v — (a AV) - (S30)v, 


where pm = pm represents the proper mass density. We see that if S3 = 0, we promptly 
obtain the desired expression for a stress—energy tensor that is appropriate to describe a 
pressureless ideal fluid 

T(a) =T(a) = pma- vv. (A114) 


Appendix C. GTG and Matter 


By virtue of GC, GTG is constructed such that gravitational effects are described by 
a pair of gauge fields, h(a) = h(a,x) and w(a) = w(a,x), defined over a flat Minkowski 
background spacetime [26], where x is the STA position vector, which is usually suppressed 
for short. 

The first of them, (a), is a position-dependent linear function mapping the vector 
argument a to vectors. The introduction of (a) ensures covariance of the equations under 
arbitrary local displacements (or an arbitrary remapping x’ = f(x)) of the matter fields in 
the background spacetime. In order to understand the physical meaning of the h field, we 
first define the covariant displacement transformation as 

M(x) 2, M(x) = M(x’), (A115) 
so that the equations A(x) = B(x) and A(x’) = B(x’) have exactly the same physical 
content. Suppose we have a vector field b(x) = V(x), where p(x) is a scalar field that is 
already covariant under displacement, i.e., p'(x) = $(x'). Now can we write b’(x) = b(x’) 
or Vo! (x) = Vyp(x')? Using the chain rule, we find 


aia (a VFO) Vx P(x’) 

a: xe WP x! 

f(a) Vap) (A116) 
a: fV’), 


where f(a) = a- = ie is a linear function of a and an arbitrary function of x, and 
f(Vw) = Vf (x) + Vx, and we call f the adjoint of f, satisfying a- f(b) = f(a) - b, or 
f(a) = d (f (b)a). It ee that Vo' (x) = f(Vyo(x')), or 


Vz =f(Vx) and PUT) HV, (A117) 


which shows us that b(x) is not covariant under displacement. In order to make objects 
such as b(x) covariant, we must introduce a position-gauge field h(a, x), which is a linear 
function of a and arbitrary function of x, so that 


nan 2, Men 2’). (A118) 


W Seg (x)) 
h(f-*(V¢"(x))) 
h(V G(x") = B(x’), 
which becomes covariant. The h(a) field plays the same role as vierbein in the tensor 
calculus approach of gauge theory of gravity [25,26]. For later use, we give the relationship 


(A119) 
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between h(a) and the metric tensor guy in GR. We define a position gauge invariant 
directional derivative as [26,28] 


Lag =a-h(V) =h(a)-V, (A120) 


where h is the adjoint of h defined by h(b) = 9,(h(c) - b), and a is an invariant vector 
(i.e. for x! = f(x), a(x) is transformed to a’(x) = a(x’)). So, h maps tangent vectors to 
tangent vectors and h maps cotangent vectors to cotangent vectors. For a given set of 
coordinates {x¥, u = 0,1,2,3}, we introduce the basis vectors 


e! = Vx", (A121) 


x 
e aad 
He gx” 


which satisfy e, -e” = On: From these vectors, we further define vectors 
8u =h (eu), gt =hler). (A122) 
These vectors satisfy the relation 
Sug = h (ey) h(e”) = ep: h1 (h(e”)) = eu: e" So). (A123) 
The metric tensor is then given by 
Suv = Qu ` gv- (A124) 
Let x(T) be a time-like curve (where T is the proper time), a mapping f : x > x’ = f(x) 
induces the transformation 


2. dx „_ dx dx 
ter e qe Y (A125) 


Comparing Equation (A125) with Equation (A120), we introduce an invariant velocity 
v = v(x(T)) as [28] 
żt=h(v), v=h1(%). (A126) 


From the known formula dx = dx"e,, the invariant normalization v? = 1 induces the 
invariant line element on a time-like curve in GR 


dt? = [Wt (dx)]? = gudx"dx". (A127) 


Another gauge field, w(a), is a position-dependent linear function mapping the vector 
argument a to bivectors. Its introduction ensures covariance of the equations of physics 
under local Lorentz rotations described by the rotor R. Under local Lorentz rotations, 
the multivector M transforms as M’ = RMR and the spinor w transforms as ’ = Ry. 
To ensure covariance of the quantities like h(V )M and h(V)w under local Lorentz rotations, 
h(V) = h(dq)a - V should be replaced by a covariant derivative D [26]. To achieve this, we 
focus attention on a- Vi’ = a - V (Ry) and write 


a-V(Rwp) = Ra- Vy + (a-VR)y. (A128) 


Clearly, the presence of the term (a - VR) renders the operator a - V non-covariant. Since 
the rotor R satisfies RR = 1, we find that 


a-VRR+Ra-VR=0, (A129) 


which implies 
a-VRR = —Ra- VR = —(a- VRR)~. (A130) 
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Hence, a - VRR is equal to minus its reverse and thus must be a bivector in the STA. We 
can therefore rewrite (A128) as 


a-V (Ry) = Ra: Vy + 5 (2a - VRR) (Ry). (A131) 


This suggests that to achieve a covariant derivative we must add a connection term to a: V 
to construct an operator 


Dap =a: VY + Talay (A132) 


which must be covariant under local Lorentz rotations. The connection Q (a) = Q(a, x) isa 
bivector valued linear function of a with an arbitrary x dependence. Under local rotations 
we expect that the operator D, will be unchanged in form, namely, 


pH =y Toa), (A133) 


Here, we have used the fact that a - V cannot change under local rotations. Nevertheless, 
the property that the covariant derivative must satisfy is 


In f 10) 
k z ee oe ae 
From these identities, it follows that O(a) transforms as 
QO! (a) = RO(a)R — 2a- VRR. (A135) 
Now, we reassemble the covariant derivative (A132) with the term h(d,) to form 
D = h(da)Da, (A136) 
and write 1 
Dy = h(d,)(a-Vpt+ 5 2(0)¥). (A137) 


Note that the vector derivative D is fully covariant, but D4 is not, since it contains the 
Q(a) field, which must transform in the same way as a - VRR under displacement and thus 
picks up a term in f (A66). Recall the definition of the position gauge invariant directional 
derivative (A120), we thus define 


a- Dy =a-h(V)p + Zola), (A138) 
where w(a) is defined by 
w(a) = O(h(a)). (A139) 


It should be pointed out that the vector a is declared to be position gauge-invariant, as stated 
below Equation (A120). Therefore, D, is related to a - D by 


Dap = 7; (a) - Dy 
1a) -A(V)p + 5w(h*(a))p (A140) 


or simply 
D, = h`! (a) - D. (A141) 


We thus also have 


Dy =a, (a h(V) p+ zew); (A142) 
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We are now ready to establish the form of the covariant derivatives of the observables 
formed from a spinor field. In general, such observables have the form 


M = 9%, (A143) 


where T is a constant multivector formed from combinations of yo, 73 and ig3. Under dis- 
placements x’ = f(x), the spinor fields are covariant: wp’ (x) = y(x"). Thus, the observable 
M inherits its transformation properties from the spinor y, 


M(x) =f, mt(x) = M(x’), (A144) 


exactly the same form in (A115). Under rotations, the spinor transforms as y’ = Ry. So it 
follows from (A143) that, under rotations, M transforms as 


M(x) È M = wT! = Ry PR = RMR. (A145) 
To achieve the fully covariant derivative for M, we first write 
LaM = Lpr ý + pT (Lap). (A146) 
This equation is a displacement covariant. We simply replace L, with a - D to obtain 
a-DM = a-Dylp+yl(a-Dy)~ 
= L(y) + (a) x (YT) (A147) 
= LaM+u(a) x M. 


Note the difference between the forms of a - D acting on M and y. 
The field strength corresponding to the w(a) gauge field is defined by 


[a-D,b- Diy = 5 Ra Ab), (A148) 
where 
R(a nb) = Law(b) — Lyw(a) + w(a) x w(b), (A149) 


a and b are constant vectors. The Ricci tensor R (a), Ricci scalar R, and Einstein tensor G (a) 
are defined, respectively, as 


R(a) = a-R(bAa), (A150) 
R = deR. (A151) 
gi = R(a) — FaR. (A152) 


The overall action integral is of the form 
1 
j= I d*x| det(h) 1(5R — Lm), (A153) 


where Lm describes the matter content and x = 87. In this paper, we adopt the covariant 
Lm = Diy3t) — mp from Dirac theory to describe the macroscopic matter, which has been 
well-studied [26,27,30] for electrons. From (A153), we obtain the following equations that 
describe the field coupled self-consistently to gravity [27]: 


torsion: DAh(a) = kh(a)-S, (A154) 


Einstein: G(a) KT (a) (A155) 
Dirac: Dypiog3 = mpyo, (A156) 
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where D ^ h(a) is the gravitational torsion which is determined by the trivector spin density 
= 1piy, x = 87, and 


T(a) = (a: Dipiys)1 (A157) 
is the matter stress—energy tensor. We can solve Equation (A154) for w(a) to obtain [27] 
' 1 1 1 
w(a) =w'(a) + aka S = —H(a)+ ha [dp A H(b)] + aKa: S, (A158) 


this defines w’ (a) as the w-function in the absence of torsion, and 


H(a) =h(V nht (a)) = —A(V) a h(a (a)), (A159) 


where ‘overdot’ notation is employed to denote the scope of a differential operator. This 
proves convenient if we employ the primed symbols to denote the torsion-free part of the 
curvature tensors, and we obtain [27] 


R(anb) = Ri(adb)+ ella Ab)-S]-S 
—5x[(a Ab) -D|-S, (A160) 
R(a) = R'(a)+ L(a -S)- S 
-ixa -(D-S), (A161) 
R = R+ g2, (A162) 
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